NATIONAL ADVISORY COMMITTEE 


FOR AERONAUTICS 


REPORT 1292 


ie Garrett Cnrpsraaon 

AiResearch Mfg. Div. 

A U G 2 1 - ■ 

LIBRARY 


INTENSITY, SCALE, AND SPECTRA OF TURBULENCE 
IN MIXING REGION OF FREE SUBSONIC JET 


By JAMES C. LAURENCE 



1956 


For sale by the Superintendent of Documents, U. S. Government Printing Office, Washington 25, D. C. Yearly subscription, $10; foreign, $11.25; 

single copy price varies according to size - Price 30 cents 


REPORT 1292 


INTENSITY, SCALE, AND SPECTRA OF TURBULENCE IN 
MIXING REGION OF FREE SUBSONIC JET 

Bv JAMES C. LAURENCE 


Lewis Flight Propulsion Laboratory 
Cleveland, Ohio 


i 


National Advisory Committee for Aeronautics 

Headquarters , 1512 II Street A 7 ir., Washington 25 , ZX C. 

Created by act of Congress approved March 3, 1915, for the supervision and direction of the scientific study 
of the problems of flight (U. S. Code, title 50, sec. 151). Its membership was increased from 12 to 15 by act 
approved March 2, 1929, and to 17 by act approved May 25, 1948. The members are appointed by the President, 
and serve as such without compensation. 


Jerome C. Hunsakeu, Sc. D., Massachusetts Institute of Technology, Chairman 
Leonard Carmichael, Ph. D., Secretary, Smithsonian Institution, Vice Chairman 


Joseph P. Adams, LL.B., Vice Chairman, Civil Aeronautics 
Board. 

Allen V. Astin, Ph. D., Director, National Bureau of Standards. 
Preston R. Bassett, M. A., Vice President, Sperry Rand Corp. 
Detlev W. Bronk, Ph. D., President, Rockefeller Institute for 
Medical Research. 

Thomas S. Combs, Vice Admiral, United States Navy, Deputy 
Chief of Naval Operations (Air). 

Frederick C. Crawford, Sc. D., Chairman of the Board, 
Thompson Products, Inc. 

James H. Doolittle, Sc. D., Vice President, Shell Oil Co. 
Clifford C. Furnas, Ph. D., Assistant Secretary of Defense 
(Research and Development) Department of Defense. 


Carl J. Pfingstag, Rear Admiral, United States Navy, Assistant 
Chief for Field Activities, Bureau of Aeronautics. 

Donald L. Putt, Lieutenant General, United States Air Force, 
Deputy Chief of Staff (Development). 

Arthur E. Raymond, Sc. D., Vice President — Engineering, 
Douglas Aircraft Co., Inc. 

Francis W. Reichelderfer, Sc. D., Chief, United States 
Weather Bureau. 

Edward V. Rickenbacker, Sc. D., Chairman of the Board, 
Eastern Air Lines, Inc. 

Louis S. Rothschild, Ph. B., Under Secretary of Commerce 
for Transportation. 

Nathan F. Twining, General, United States Air Force, Chief 
of Staff. 


Hugh L. Dryden, Ph. D., Director 
John W. Crowley, Jr., B. S., Associate Director for Research 


John F. Victory, LL. D., Executive Secretary 
Edward H. Chamberlin, Executive Officer 


Henry J. E. Reid, D. Eng., Director, Langley Aeronautical Laboratory, Langley Field, Va. 
Smith J. DeFrance, D. Eng., Director, Ames Aeronautical Laboratory, Moffett Field, Calif. 
Edward R. Sharp, Sc. D., Director, Lewis Flight Propulsion Laboratory, Cleveland, Ohio 
Walter C. Williams, B. S., Chief, High-Speed Flight Station, Edwards, Calif. 


ii 


REPORT 1292 


INTENSITY, SCALE, AND SPECTRA OF TURBULENCE IN MIXING 
REGION OF FREE SUBSONIC JET 1 


By James C. Laurence 


SUMMARY 

The intensity of turbulence, the longitudinal and lateral corre- 
lation coefficients, and the spectra of turbulence in a S. 5 -inch- 
diameter free jet were measured with hot-wire anemometers at 
exit Mach numbers from 0.2 to 0.7 and Reynolds numbers from 
192,000 to 725,000. 

The results of these measurements show the following: ( 1 ) 
Near the nozzle {distances less than 4 or 5 jet diam downstream 
of the nozzle) the intensity of turbulence , expressed as percent 
of core velocity, is a maximum at a distance of approximately 
1 jet radius from the centerline and decreases slightly with 
increasing Mach and for Reynolds number. At distances greater 
than 8 jet diameters downstream of the nozzle, however, the 
maximum intensity moves out and decreases in magnitude until 
the turbulence-intensity profiles are quite flat and approaching 
similarity. {2) The lateral and longitudinal scales of turbu- 
lence are nearly independent of Mach and/or Reynolds number 
and in the mixing zone near the jet vary proportionally with 
distance from the jet nozzle. {3) Farther downstream of the 
jet the longitudinal scale reaches a maximum and then decreases 
approximately linearly with distance. {4) Near the nozzle the 
lateral scale is much smaller than the longitudinal and does not 
vary with distance from the centerline, while the longitudinal 
scale is a maximum at a distance from the centerline of about 
0.7 to 0.8 of the jet radius. {5) Farther downstream this maxi- 
mum moves out from the centerline. {6) A statistical analysis 
of the correlograms and spectra yields a “ scale ” which , although 
different in magnitude from the conventional , varies similarly 
to the ordinary scale and is easier to evaluate. 

INTRODUCTION 

Recent analyses by Light-hill (ref. 1) and by others (refs. 2 
and 3) show that turbulence may be chiefly responsible for 
the noise from high-speed jets. However, a complete analy- 
sis relating turbulence to noise has not yet been made. 
Consequently, it is not known where the noise originates in 
the jet or which of the turbulence parameters, intensity, 
scale, or spectral distribution, is best suited to noise studies. 
Lighthill has suggested that each turbulent eddy might be 
considered as a single concentrated sound radiator. If 
Ligh thill’s suggestion is valid, it will be necessary to know 
the turbulence characteristics throughout the jet. The 
intensity of sound radiation will surely he related to the 
intensity of turbulence. Furthermore, the spectral distribu- 


tion of turbulent energy might reasonably be expected to 
give information about the spectral distribution of sound 
energy. 

Three regions of the jet need to be investigated (see fig. 1): 



(1) the central core, where the velocity profiles are flat and 
the intensity of the turbulent fluctuations is low; (2) the 
mixing region, which lies between the central core and the 
undisturbed air of the surroundings, where the velocity 
gradients are large and the intensity of turbulence is high; 
and (3) the region farther downstream where the central core 
and the mixing region blend into a completely turbulent air- 
stream. The purpose of experiments in these regions is to 
measure the statistical parameters of the turbulence in the 
mixing zone and in the completely turbulent region of the jet. 

Several investigators have reported the results of work in 
these regions (refs. 4 to 8). Reference 4 has only a limited 
amount of hot-wire data in the mixing zone, along with 
velocity-profile and turbulence-level measurements near the 
axis of the jet. Corrsin has made or has collaborated in 
making the most complete measurements in jets to be found 
in the literature (refs. 5 to 7). These investigations are of 
great interest and of much use, but the velocities investigated 
are small (<100 ft/sec). Hot-wire measurements in a 
unique form of a half-jet in which a solid wall replaces one 
side of the free mixing region are reported in reference 8. 

In order to obtain data in speed ranges that are of more 
interest in jet noise research, a systematic investigation of 
the turbulent mixing region of a free jet was undertaken as 
part of the jet noise program at the NACA Lewis laboratory. 


i Supersedes NACA Technical Notes 3561, “Intensity, Scale, and Spectra of Turbulence in Mixing Region of Free Subsonic Jet/’ by James C. Laurence, 1955; and 3576, “Further Measure- 
ments of Intensity, Scale, and Spectra of Turbulence in a Subsonic Jet,” by James C. Laurence and Truman M. Stickney, 1956. 


1 


2 


REPORT 1292 — NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 


Principal interest was placed on measurements with the hot- 
wire anemometer of the statistical parameters of turbulence 
in the exit Mach number range from 0.2 to 0.7 and Reynolds 
numbers (based on jet radius) from 192,000 to 725,000. 
Since it was impossible to vary Mach numbers and Reynolds 
numbers independently, their combined effects are called 
Mach and/or Reynolds number effects herein. 

The analysis of the results of these turbulence studies 
follows the statistical approach first introduced by Taylor 
(ref. 9). Dryden reviews this method in reference 10 and 
shows the success of this type of approach in dealing with 
isotropic turbulence. Similar approaches to the analysis of 
the results of studies of jet turbulence have been made also. 
But the results have not been very successful, especially in 
the mixing region of the jet. Hence, a statistical analysis is 
presented in this report which uses a different value of the 
“scale” or eddy size based on the energy spectrum of the 
turbulence. 


INSTRUMENTATION AND TEST FACILITIES 

HOT-WIRE ANEMOMETERS 

The hot-wire anemometers used in this test were the 
(‘distant -temperature anemometers described in reference 1 1 . 
Constant-temperature anemometers were chosen instead of 
constant-current anemometers because of certain advantages: 
(1) The wire current is controlled by an electronic servosys- 
tem that protects against accidental wire burnout when the 
flow is suddenly reduced. (2) Since the bridge balance is 
maintained during flow changes, there are no compensation 
controls to set and no testing to determine whether the com- 
pensation is effective. (3) Finally, this instrument can fol- 
low large fluctuations in flow without appreciable error. 

The hot-wire anemometer and the auxiliary equipment- 
have been improved in several respects (see refs. 1 1 and 12). 
The direct-current amplifier has an improved frequency re- 
sponse, and its equivalent input noise has been decreased to 



(a) Longitudinal correlation of velocity components in z-direction. 

(b) Lateral correlation of velocity com- (c) Intensities and autocorrelations of (d) Velocity components in ^-direction, 

ponents in ^-direction. velocity in ^-direction. 


Figure 2. -Hot-wire-anemometer probes. 
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a value comparable with that of other hot-wire instruments 
(less than 100 microvolts). Frequency-response measure- 
ments made according to the methods described in reference 
1 1 show an over-all response (amplifier, bridge, wire, and 
cables) of about 100 kiloc 3 ~cles. 

HOT-WIRE PROBES 

The hot-wire probes (fig. 2) were single- and double-wire 
probes. The double-wire probes were arranged in parallel 
and X-array. The parallel arrays (figs. 2(a) and (b)) were 
used to measure the longitudinal and lateral u velocity cor- 
relation coefficients, respectively. For the longitudinal cor- 
relations one wire was held fixed and the other was moved 
downstream. This wire was displaced not more than 0.005 
inch towards the centerline of the jet to avoid interference 
from the flow over the first wire. A similar arrangement was 
used for the lateral correlations. Practical limitations on the 
separability of the wires limited the number of points for 
which measurements could be obtained. The measured cor- 
relations were assumed to be those at the position of the 
fixed wire. The X-array (fig. 2(c)) was employed to meas- 
ure the ^-component of the velocity fluctuations. The 
single-wire probe (fig. 2(d)) measured the u velocity in- 
tensities and all autocorrelations. 

The wire material used in the tests was the 0.0002-inch 
tungsten wire evaluated in reference 11. The mounting 
procedure is an adaptation of the plating-soldering technique 
that was first used by the National Bureau of Standards 
(ref. 13). 

SPECTRUM ANALYZERS 

Two spectrum analyzers were used in obtaining the spectral 
density data. One of these has a constant narrow band pass 
(5 cps) with a frequency response from 20 to 16, 000 cps; 
the other is a constant-percentage band-pass instrument 
with about the same frequency response. The filter char- 
acteristics of the two instruments are given in figure 3. The 
two analyzers give comparable spectral density curves within 
the general accuracy of the experiments. 



TEST FACILITIES 


TAPE RECORDER 

A dual-channel tape recorder was used to record the signals 
from a single hot-wire-anemometer probe located in the 
airstream. These recordings were later used to obtain the 
autocorrelation coefficients. 

ACTUATORS 

The lateral positions of the probe were controlled by screw- 
type actuators. The longitudinal positions, however, were 
set by adjustments of a lathe carriage. The probes were 
alined in the stream with reference to the centerline of the 
jet. This aerodynamic centerline was determined from 
total-pressure measurements in the jet. A plug was ma- 
chined to fit the exit nozzle of the jet, and a steel wire was 
stretched along the centerline to a point on the centerline 
far downstream in the jet (approx. 12 ft). A telescope 
mounted at this point helped in relocating the probes after 
a probe failure without restringing the steel wire. 


The subsonic jet and associated plenum chamber used in 
these experiments are shown schematically in figure 4. 
The jet has a 3.5-inch-diameter nozzle at the end of a bell- 

Compressed-air 

supply 



Figure 4. — Schematic diagram of subsonic jet, plenum chamber, 
filters, and location of pressure- and temperature-measuring in- 
struments. 
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mouth which is attached to a large plenum chamber. The 
turbulence intensity in the plenum chamber varied from 4 
to 5 percent over the range of the test conditions. The 
whole is connected to a compressed-air supply. A filter is 
included as an integral part of the plenum chamber. It 
successfully removes a large part of the dust and pipe scale 
sometimes found in similar airstreams. The flow is adjusted 
by means of electrically controlled valves throughout the 
subsonic range of Mach numbers. The instruments for 
measuring total pressure and temperature are placed in the 
plenum chamber and read remotely. Figure 1 shows the 
approximate location of the testing points in the jet, and 
table I gives the exact coordinates. 


TABLE I. — LOCATION OF TEST POINTS IN SUBSONIC JET 


Distance 
from jet 
nozzle, 
x/r 



Distance from jet centerline. 

y/r 



1. 14 

1.206 

1. 154 

1. 103 

1.051 

1.00 

0. 949 

0.846 

0. 794 

0.743 

2.29 


1.308 

1.206 

1. 103 

1. 00 

.897 

. 794 

.692 



4.58 


1.411 

1. 200 


1. 00 

.794 

.589 

.383 


7. 0 


1.686 

1. .500 

1*343 

1.00 

. 657 

.314 

-.029 


8 


2.000 

1.714 

1.280 

1. on 

.714 

.280 

0 


12 


2. (K) 

1.714 

1.280 

1.00 

.714 

. 280 

0 



16 


2. 00 

1.714 

1.286 

1.00 

.714 

. 286 

0 



24 

5’ 143 

4.400 

3. 772 

3.088 

2. 400 

1.715 

1.028 

.343 

0 

32 

5. 143 

4. 400 

3. 772 

3. 088 

2. 400 

1.715 

1.028 

.343 

! 0 

40 

5. 143 

4. 400 

3. 772 

3.088 

2. 400 

1.715 

1.028 

.343 

1 0 


EXPERIMENTAL PROCEDURE 
INTENSITY MEASUREMENTS 

A single hot wire was used to measure the ^-component 
of the turbulent intensity. An X-wire probe was used to 
measure the ^-component. The hot-wire signal (or the 
instantaneous difference between the two X -wires) was meas- 
ured directly by a root-mean-square voltmeter. Reference 
1 1 describes the root-mean-square voltmeter (average- 
square computer) which was used in an improved form for 
the first half of the intensity surveys made. A true rms 
voltmeter was used for the remainder of the data. It gives 
a precise wide-band squaring action that is unaffected by 
normal ambient temperature variations and, in general, is 
more satisfactory than the other instrument. 

The turbulence intensity was calculated by methods 
outlined in references 11, 13, and 14. 

CORRELATION MEASUREMENTS 

Velocity correlations in space and in time are of interest 
in the noise problem. Two hot wires can be used to obtain 
the lateral and longitudinal correlations. The wires are 
mounted on telescoping tubes that can be displaced either 
laterally or longitudinally with respect to each other by 
actuators. During the actual measurements, the fluctuating 
components of the two hot-wire signals, sayjj and e 2 , are 
led to a correlator that measures e l e 2) e?, and ef. 

If a linear relation between the velocity fluctuation and 
the hot-wire voltage is assumed, that is, 


‘K 


k\k 2 


€ 1*2 


V ^1 € 1 V ^14 


e l 6 J 

\ u \ v 


( 3 ) 

( 4 ) 


(All symbols are defined in appendix A.) 

The displacement of the wires in the lateral direction 
gives the lateral correlation coefficient; in the longitudinal, 
the longitudinal correlation coefficient. 

An autocorrelation coefficient, that is, the correlation 
between two segments of the same signal separated in time, 
can be obtained if some method is devised to obtain from the 
signal of a single hot wire another signal delayed in time 
with respect to the first. An instrument that provides time 
delays of this nature (described in refs. 16 and 17) was 
designed and built at the NACA Lewis laboratory (ref. 18). 



Figure 5. — Block diagram of correlation computer. 


In general, this correlator (see fig. 5) works as follows: 
The signal from a single hot wire is recorded simultaneously 
on two channels of a dual-channel tape recorder. A special 
playback instrument uses two pickup heads, one of which 
can be displaced with respect to the other by means of a 
micrometer screw. When the two heads are in such a posi- 
tion that they are picking up identical signals, the correla- 
tion between these signals is necessarily unity. But if the 
movable head is displaced so that it is picking up a signal 
recorded a short time earlier or later, the two signals are, 
of course, different. These two signals are the two voltages 
d and €2 considered in equation (1). The correlation 
coefficient C R is then the autocorrelation coefficient 

The value of the time interval A/, is obtained from the 
actual displacement of the head and the constant tape speed: 

( 5 ) 


U 2 =k 2 e 2 J 


as is done in reference 
signals, 


15, there results, upon combining the 


%= 


U\U 2 


( 2 ) 


It is rather obvious that there should be a relation between 
the autocorrelation and the longitudinal correlation coeffici- 
ents involving the stream velocity U and the head separation 
(ref. 19). 

An extension of this method can be made to measure the 
autocorrelation of the ^-components of the turbulence. In 
this method the difference signal from an X-wire probe 
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located in the stream is recorded simultaneously on both 
channels of the recorder and played back in the usual way. 

The double correlator of reference 1 1 was used to obtain 
correlations of the hot-wire signals for the first half of the 
program. The ratios were calculated numerically, because 
the ratio meters that are sometimes used are unreliable. 
Inaccuracy results from the use of an average-square com- 
puter, which depends upon the operating characteristics of a 
vacuum tube to obtain the square of the sum and of the 
difference of two voltages in order to apply the quarter- 
square principle of multiplication. 

After about half the data were obtained, an electronic 
multiplier replaced the double correlator. This multiplier 
(described in ref. 18) does not use the quarter-square principle 
of multiplication. The multiplier operates from d.c. to 50 
kilocycles (3-db point) and is useful in obtaining correlation 
coefficients. Two channels are available to give the average 
product of two input voltages as well as their average squares. 
The outputs can be read on any higli-impedance d-c volt- 
meter. By properly adjusting the level of the inputs, the 
correlation coefficients can be read directly. 


SPECTRUM OF TURBULENCE 


If the hot-wire signal is analyzed by means of a series of 
band-pass filters, a spectrum results. The relative energy 
spectrum is more fully discussed in references 10 and 20. 
Most of the spectra presented here were recorded on a paper 
strip chart by an analyzer of constant percentage (1/10- 
decade) band width. The switching from one filter to the 
next is done automatically and in synchronization with the 
movement of the chart. The spectrum recordings are made 
rapidly (about 1 per min) and with little attention by the 
operator. 

The chart recordings are in decibels above any practical 
reference voltage, and at one point on the chart the total 
voltage level for the entire range of frequencies is recorded. 
The conversion to spectral density is as follows: By definition, 


Therefore, 


N b =20 log^- & 

tref 


iV,„,=20 log 5= 


tref 


antilog § 


*ref 


10 


antilog ^ 


*ref 


10 


( 6 ) 


(7) 


Hence, the spectral density function 


iV s 


antilog 

y</>4 


antilog —j 1 


( 8 ) 


where B is the band-pass width of the filters. 

The resulting display of the spectral density function with 
frequency shows the content of the signal in each individual 
frequency band. This method of analysis is important foi* 
two reasons: (1) It tells the experimenter the frequency bands 


in which most of the energy is concentrated, and (2) the 
spectrum can be used to obtain the average eddy size. 

SCALE OF TURBULENCE 

The scales of turbulence (ref. 21) are defined in many 


applications as follows: 


Longitudinal: 



(9a) 

Lateral : 



(9b) 

r» oo 

(9c) 



These scales are generally regarded as the physical dimen- 
sions in the x - (longitudinal), y - (radial), and 2 - (tangential) 
directions of the average eddy in the flow. The three scales, 
as defined in equations (9), are all important in a study of jet- 
noise. However, variations in the tangential direction are 
not included in this report. Reference 21 shows that the 
correlation coefficient is related to t he spectral density as 
follows : 

J(f)cOS~^6f ( 10 ) 

and, inversely, 

J {f) =M»^ x cos nr dx (11) 

Thus, there are two possible ways to obtain the longitudinal 
correlation coefficient. One is to measure it directly; the* 
other is to obtain it from the turbulence spectrum and 
equation (10). If the longitudinal correlation coefficient 
is assumed to be exponential in form, that is, 

%=e~ xlJ:x (12) 

then the transform (eq. (11)) can be readily worked out as 

(ref. 10) 



Analogous relations are available for the transverse correla- 
tions. 
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The methods just outlined have been used to great ad- 
vantage in the analysis of turbulence induced by grids in 
wind tunnels (ref. 21). Similar application to the flow in 
subsonic jets does not prove to be so successful, however, 
particularly in the mixing region close to the nozzle of the 
jet, since the definitions do not apply to shear flows. 

In particular, the correlograms that resulted from the two- 
wire measurements and the autocorrelation measurements in 
the jet are not fitted very well by a single exponential. A 
typical correlogram is shown in figure 6(a). It is apparent 
that the shape of this curve is different from a pure ex- 
ponential. In order to find a curve of best fit for this data, a 
function was built up which is the algebraic sum of two ex- 
ponentials, as shown in figure 6(b). The correlation fitting 
function is then of the following form (see appendix B): 



l-^ 

(14a) 

% 

~ u*t 

1-0 

(14b) 


In these equations the constants /h A x , and w* have special- 
ized meanings. The A x is a characteristic eddy size that is 
different in magnitude from the longitudinal scale of the tur- 
bulence ^ r , and co* is related to this eddy size by co*=t//A J . 
The parameter @ is related to the cut-off frequencies (upper 
and lower) of the spectral density curves by the relation 



(a) Typical autocorrelogram. 

(b) Method used to approximate aut.ocorrelograms. 

Figure 6. — Autocorrelograms; /3 = 0.2; co* = 6283 radians per second. 


fu (1+0) 2 

The cut-off frequency is defined in the usual electronic sense 
as the intersection of the asymptotes of the spectral density 
curve with a horizontal line through the maximum point of 
the spectral density curve. The relation between the distance 
jc and the time t in these experiments is the velocity of the 
stream; that is (ref. 20), 

x=U l0C t (15) 

Equations (14) can be adjusted by means of the quantities 
A x , 0, and co*to fit the correlograms obtained experimentally. 

It is now obvious that the conventional definition of the 
scales of t lie turbulence (eqs. (9)) is no longer valid, since the 
integral is always zero if has the form of equations (14). 
The parameter however, is useful in obtaining a measure of 
the eddy size, which, after all, is what is meant by the scale 
of the turbulence. In equation (14a), if 7^ z = 0, 


_j o (0-n 

In 0 


(16) 


where j^o is the value of the wire separation when is first 
equal to zero. This definition of scale is related to two quan- 
tities that are quite easy to evaluate from the experimental 
correlograms and spectra: (1) the first crossing of the abscissa 
by the correlation curve, and (2) the parameter f$, which is 
related to the ratio of t lie low- to the high-frequency cut-off 
of the spectral density curves. 



Figure 7. — Typical spectral density curve. 
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The transform of equation (14a) can be obtained with the 
following result: 


tv 1 T*j3 

7(/)=4fT 






(17) 


where /3 is the same as in equations (14), and 



(18) 


Figure 7 siiows a typical spectral density curve obtained in 
the completely turbulent region of a subsonic jet. It repre- 
sents data obtained at the same time and at the same location 
in the same jet as the correlogram of figure 6(a). A spectrum 
fitting function of the type of equation (17) was calculated by 
the methods described in appendix B from the correlation 
curve. These calculated points are shown for comparison 
with the experimentally obtained spectrum. 

ACCURACY OF MEASUREMENTS 

The accuracy of the hot-wire measurements made with the 
constant-temperature instrument is discussed in references 1 1 
and 12. It is sufficient to say here that these measurements 


are made as accurately as the usual hot-wire measurements. 
In measuring quantities that vary randomly with time, a 
certain amount of discretion must be used in obtaining aver- 
aged readings of meter pointers that are jumping about. 
Reference 12 includes an evaluation of the NACA Lewis 
laboratory hot-wire equipment and experimental techniques 
and compares them with other systems. 

No corrections have been made to any of the hot-wire re- 
sults for the finite length of the wire. In these measurements, 
the scales of turbulence were much greater than the length 
of the hot wires (0.080 and 0.040 in.). 

Some of the intensities measured were large compared with 
the local mean flow. In the constant-temperature method of 
anemometry, these large intensities can be evaluated (see ap- 
pendix F, ref. 11). In the usual method and in the method 
used herein, however, the root mean square of the voltage 
fluctuations is used in place of their instantaneous values. 
This approximation is valid only for small values of the fluctu- 
ations. Nevertheless, vibration studies (e. g., ref. 5) have 
shown that the hot-wire anemometer itself is capable of fol- 
lowing fluctuations in a single direction faithfully up to 60 
to 70 percent. 

When the fluctuations in U become large, errors will be 
introduced because of the interaction of the components of 
the turbulence; that is, the effect of u , v , and w on U and on 



Distance from jet centerline, y/r 


(a) Distance from jet nozzle x/r, (b) Distance from jet nozzle x/r, (c) Distance from jet nozzle x/r, (d) Distance from jet nozzle x/r, 
1.14. 2.29. 4.58. 7.60. 

(e) Distance from jet nozzle x/r, 7.60, 8, and 12. Exit Mach number, (f) Distance from jet nozzle x/r, 16 to 40. Exit Macli number, 0.3; 
0.3; Reynolds number, 300,000. Reynolds number, 300,000. 

Figure 8. — Intensity of turbulence in percent of core velocity at various exit Mach and/or Reynolds numbers. 

398434—57 2 
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each other. These errors are considered in reference 12, 
where actual corrections are worked out. No such correc- 
tions have been applied to the results shown in this report. 

Some of the correlations measured with the average-square 
computer and the sum and difference circuits (quarter-square 
multiplication) show values greater than unity. This is, of 
course, impossible, and such readings must be the result of 
inaccuracies in the multiplication. In most of these in- 
stances enough points are available to fair a reasonable curve 
without considering these points (7^^>l). 

The largest sources of error in these experiments were (1) 
the inaccuracy in repeating setups after probes were changed, 
(2) the fluctuations in pressure in the air supply, which 
caused fluctuations in the mean flow level, and (3) the 
presence in the flow of periodic fluctuations due to sound or 
pressure waves and floor vibrations. It is believed that the 
measurements reported herein for turbulence intensities 
loss than 30 percent are not in error by more than ±10 
percent. 

RESULTS AND DISCUSSION 

The results of these hot-wire measurements are presented 
in a series of graphs and tables. In general, the points for 
which the measurements are reported are shown in figure 1 
and are tabulated in table I. However, for the intensity 
surveys many more points were used. For some of the 
correlation measurements, physical limitations in the amount 
of wire separation made it impossible to obtain results at 
all the points indicated in figure 1. 

The order in which the results are presented has no special 
significance. The surveys are given for only one-half of the 
jet, because the jet was found to be essentially axisymmetric. 
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Distance from jet centerline y/r- 0, percent of mixing-zone radius 

(a) Distances from jet nozzle less than 12 jet radii. 

(b) Distances from jet nozzle greater than 24 jet radii. 

Figure 9.— Similarity of turbulence-intensity profiles. Exit Mach 
number, 0.3; Reynolds nurr.ber, 300,000. 


INTENSITY OF TURHULENCE 

The intensity measurements were made at values of the 
nondimensional radius yjr distributed across the mixing 
zone for several positions downstream of the jet nozzle. 
The results of these measurements are shown in figures 8 to 
10. In figures 8(a), (b), (c), and (d), the turbulence in- 
tensity in percent of central-core velocity is presented for 
x/r values of 1.14, 2.29, 4.58, and 7.60, respectively, for four 
Mach numbers (0.2, 0.3, 0.5, and 0.7). In figures 8(e) and 
(f), the intensity is shown only for a Mach number of 0.3 at 
x/r values of 7.60, 8, 12, 16, 24, 32, and 40. 

One point of particular interest is shown by figure 9. 
Figure 9(a) shows that the turbulence-intensity profiles near 
the jet nozzle (where the potential core is still present) are 
similar in shape and little different in actual magnitude of 
the intensity of turbulence. As the distance from the nozzle 
increases, however, this similarity disappears only to de- 
velop again in a different form far downstream. Figure 9(b) 
shows the reappearance of a degree of similarity, especially 
at points removed from the centerline of the jet. Here the 
turbulent velocity profiles are becoming quite flat for a 
large part of the mixing zone, the maximum value of the 
intensity moving out from the jet centerline. Figures 10(a) 
to (d) show the intensity of turbulence in percent of local 
mean velocity for the same Mach numbers and the same 
x/r values as figures 8(a) to (d). 



1.14. 4.58. 


(b) Distance from jet nozzle x/r, (d) Distance from jet nozzle x/r, 
2.29. 7.60. 

Figure 10. — Intensity of turbulance in percent of local velocity at 
various exit Mach and/or Reynolds mumbers. 
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The mean velocity profiles (obtained from total-pressure surveys) are shown in figure 11. The local mean velocity 
in terms of central-core velocity is plotted against the nondimensional radius y/r for the entire range of Mach numbers 
and x/r values up to 40. 





(a) Distance from jet nozzle x/r , 1.14 
(c) Distance from jet nozzle x/r, 4.58. 
(e) Distance from jet nozzle x/r, 8. 

(g) Distance from jet nozzle x/r, 16. 
(i) Distance from jet nozzle x/r , 32. 


(b) Distance from jet nozzle x/r, 2.29 
(d) Distance from jet nozzle x/r, 7.60. 
(f) Distance from jet nozzle x/r, 12. 
(h) Distance from jet nozzle x/r, 24. 
(j) Distance from jet nozzle x/r, 40. 


Figure 11. — Mean velocity profiles. 
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Figure 13. — Effect of potential-core Mach number on boundaries of 
3.5-inch-diameter jet. 



Figure 14. — Variation of intensity of turbulence with airspeed for 
circular nozzle. Distance from nozzle x/r, 4.58. 


The moan velocity profiles are also similar, as is shown in 
figure 12, where the ratio of the local mean velocity to the 
maximum mean velocity is plotted against a nondimensional 
mixing-zone width. This width is defined as the ratio of 
the distance from the jet centerline to the point where the 
local velocity is half the maximum velocity. The figure 
shows very well the degree of similarity for all distances 
from the jet nozzle as well as for the different Mach numbers. 
Also shown is a comparison with the theoretical momentum- 
transfer results as given by reference 22. 

The result of these mean velocity studies is shown in 
figure 13, where the approximate jet boundaries are shown 
for the different Mach numbers. The positions of the zero 
velocity, the one-half maximum, and the core velocity are 
shown. 

A comparison of figures 8 and 1 1 shows that the maximum 
intensity occurs near the points of inflection of the velocity 
profiles at approximately 1 jet radius from the centerline 
for distances from the jet nozzle up to about 10 jet radii. 
As the distance from the nozzle increases, however, this posi- 
tion of the maximum intensity begins to move out from the 
jet centerline, and similarity of the profiles tends to develop 
far downstream from the nozzle. The curves show that the 
line of maximum shear moves out from the centerline as the 
distance downstream of the jet nozzle increases. 

The observed turbulence levels decreased with increasing 
Mach and/or Reynolds number throughout the jet. Figure 
14 shows that the level of turbulence in the central core and 
at the position of maximum intensity decreased with increas- 
ing exit Macli and/or Reynolds number. It is not clear, 
however, that the turbulence which is generated in the jet 
mixing process would change with jet Mach number if there 
were no initial turbulence in the jet. The decrease in inten- 
sity shown for the core may possibly explain the decrease 
noted in the mixing zone. If this explanation is accepted, 
the results of the intensity surveys are in a fair agreement 
with those reported in references 5 and 8. 

The turbulence intensity of the ^-component of the veloc- 
ity fluctuations was measured at four points in the jet. 
These measurements, along with the corresponding u inten- 
sities, are included in table II. This table shows that the 
u and v intensities are not equal but may differ by about 
as much as 17 percent. 


TABLE II.— COMPARISON OF u AND v INTENSITY 
MEASUREMENTS 


x/r 

y/r 

V " 2 

^ loc 

V l ' 2 

Ulo* 

Difference a , 
percent 

6 

1.00 

0. 241 

0.242 

0.4 

6 

0 

. 101 

. 118 

10.8 

8 

1.00 

. 266 

.241 

-10.4 

8 

0 

. 148 

. 134 

-10.4 


11 Difference is in percent of smaller intensity. 
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(a) Distance from jet centerline y/r, 1.065. 

(b) Distance from jet centerline y/r, 1.00. 

(c) Distance from jet centerline y/r, 0.935. 

Figure 15. — lateral velocity correlations. Distance from jet nozzle 
x/r, 1.14; exit Mach number, 0.3: Reynolds number, 300,000. 

LATERAL CORRELATION COEFFICIENTS 

For each distance from the jet nozzle x/r, three different 
distances from the jet centerline y/r were used in the meas- 
urement of lateral correlation coefficients for all the Mach 
numbers investigated At positions downstream greate • than 
approximately 1 jet diameter, the wires on the probes could 
not be separated far enough to reach zero correlation. Fig- 
ures 15 and 16 are typical of the results for the lateral cor- 
relation coefficients. Figure 15 shows the variation of the 
lateral correlation for x/r of 1.14 and for three y/r values. 
Figure 16 shows the lateral correlation for y/r of 1.00 at 
several x/r values. Also shown on the curves of figure 16 
are exponential curves that have been fitted to the data. 
The agreement is good for x/r values of 7.60 and 4.58 but 
only fair for x/r of 2.29 and 1.14. In some cases, however, 
this agreement justifies the use of equation (12) to obtain 
equation (13) for the spectral density function. 





0 .04 .08 .12 .16 .20 .24 .28 .32 .36 .40 .44 .48 .52 


Wire separation, 77, in. 

Figure 16. — Lateral velocity correlations. Distance from jet center- 
line y/r, 1.00; exit Macli number, 0.3; Reynolds number, 300,000. 
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Distance from jet nozzle, x/r 

Figure 17. — Variation of lateral scale with distance from jet nozzle. 
Distance from jet centerline y/r, 1.00; exit Mach number, 0.3; 
Reynolds number, 300,000. 
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Figure 18.— Effect of exit Mach and/or Reynolds number on lateral velocity correlations. Distance from 
jet nozzle x/r , 0.57; distance from jet centerline y/r, 1.00. 
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i i i J , 

192 300 500 725x10° 

Reynolds number, Re 

Figure 10. Effect of exit Mach and/or Reynolds number on lateral 
scale of turbulence. Distance from jet centerline y/r, 1.00. 

LATERAL SCALE OF TURBULENCE 

The lateral correlation data of figures 15 and 16 for 
y/r= 1.00 are compared in figure 17, where the lateral scale 
of turbulence is shown as a function of distance downstream 
of the nozzle. An exponential curve was faired through the 
data points of the correlation curves (fig. 16), and a repre- 
sentative value of the scale was obtained. The results show 
an approximately proportional increase of lateral scale with 
distance from the jet nozzle. The least-square line gives 
the relation 

oC,=0.036z+0.043 (19) 

which compares with 

-C y =0.028z 


given in reference 8. Little, if any, variation of the lateral 
scale was found with distance from the jet centerline. 

MACH AND OR REYNOLDS NUMBER EFFECTS 

In figure 18 are shown the effects of exit Mach and/or 
Reynolds number on the lateral correlation coefficient as 
a function of wire separation for Mach numbers of 0.2, 0.3, 
0.5, and 0.7. These curves show no variation of lateral cor- 
relation coefficient with Mach and/or Reynolds number and 
show clearly that there is little or no effect of Mach and/or 
Reynolds number on the lateral scale. Figure 19 summa- 
rizes the data of figure 18 and similar data for x/r of 0.57 
and 1.14 at a fixed y/r of 1.00. 

LONGITUDINAL CORRELATION COEFFICIENTS 

The longitudinal correlations measured with two wires are 
given in figures 20 and 21. Figure 20 shows the variation 
of the longitudinal correlation coefficient with distance 
across the mixing zone for x/r of 1.14. At the points near the 
core of the jet (y/r< 1), strong periodic signals of approxi- 
mately 8000 cps were encountered. Their presence causes 
the apparent scatter in the data of figures 20(g) to (i). The 
source of these disturbances is unknown, but the wavelength 
of the disturbance is approximately half the diameter of the 
jet nozzle. Perhaps transverse oscillations are responsible. 
Figure 21 shows the variation of the longitudinal correlation 
coefficient with distances from the jet nozzle of 1.14, 2.29, 
and 4.58. 


( 20 ) 


Longitudinal correlation coefficient, # x Longidinal correlation coefficient, 
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Wire separation, in. 

(a) Distance from jet centerline y/r, 1.20 (b) Distance from jet centerline y/r, 1.154. (c) Distance from jet centerline y/r, 1.103. 

(d) Distance from jet centerline y/r , 1.051. (e) Distance from jet centerline y/r, 1.00. (f) Distance from jet centerline y/r, 0.949. 

(g) Distance from jet centerline y/r , 0.897. (h) Distance from jet centerline y/r, 0.846. (i) Distance from jet centerline y/r, 0.794. 

Figure 20. — Longitudinal velocity correlations. Distance from jet nozzle x/r 1.14; exit Mach number, 0.3; Reynolds number, 300, COO. 





Inspection of figures 20 and 21 shows that the longitudinal 
correlation is affected by varying x/r and y/r. This effect 
will become more apparent when the longitudinal scale is 
evaluated from the area under the correlation curves; the 
variation is discussed in the section Longitudinal Scale 
of Turbulence. 

AUTOCORRELATIONS 

The autocorrelation functions are presented in figures 22 
and 23. Figure 22 shows the autocorrelation coefficient as 
a function of delay time in milliseconds for a fixed x/r value of 
1.14 and for y/r values distributed as shown in figure 1. 
Figure 23 shows curves of autocorrelation coefficient against 
delay time for a fixed y/r of about 1.00 for several x/r values. 
Again, the effect of varying x/r and y/r will become apparent 
when the scale of turbulence is discussed. In figures 22 and 
23 the influence of the periodic disturbances is exaggerated 
because of the size of the time-delay intervals. The 8000-cps 
sound wave is especially prominent, as well as another type 
of extraneous signal that resulted from floor vibrations due 
to heavy machinery in the building. For example, figures 
22(g) and (h) show both types of interference. These 
wiggles are not present in figure 23, where the distribution 
with distance downstream of the nozzle of the jet is shown. 

As has been pointed out, the autocorrelations can be con- 
verted to longitudinal correlations by certain relations of 
the mean flow. Figure 24 shows a comparison of the two 
measurements of the longitudinal correlation coefficient. 
The agreement is quite good, and the curves shown are 
typical of the measurements. The two-wire correlations 


Figure 21. — Longitudinal velocity correlations. Distance from jet 
centerline y/r, 1.00; exit Mach number, 0.3; Reynolds number, 
300,000. 


Long itudinal correlation coefficient, st x __ Autocorrelation coefficient, 
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(a) Distance from jet centerline y/r, 1.206. (d) Distance from jet centerline y/r, 1.051. (g) Distance from jot centerline y/r, 0.897. 

(b) Distance from jet centerline y/r, 1.154. (e) Distance from jet centerline y/r, 1.00. (h) Distance from jet centerline y/r, 0.846. 

(c) Distance from jet centerline y/r, 1.103. (f) Distance from jet centerline y/r , 0.949. (i) Distance from jet centerline y/r, 0.794. 

Figure 22. — Autocorrelations. Distance from jet nozzle x/r, 1.14; exit Mach number, 0.3; Reynolds number, 300,000. 




wore measured with the quarter-square method of multipli- 
cation using the average-square computer. The one -wire 
results were measured with the analog multiplier. If both 
measurements bad been made with the analog multiplier, 
the agreement would probably be even better. 

SPECTRA OF TURBULENCE 

The spectra of the ^-component of the turbulent velocity 
are shown in figures 25 and 26. Figure 25 shows the spectra 
for a fixed x/r value of 1.14 and a series of y/r values dis- 
tributed as shown in figure 1. On the other hand, figure 26 
presents the spectra for a fixed y/r value of approximately 
1.00 and a series of x/r values. The change in the spectra 
with distance from the nozzle of the jet is primarily a shift 
of energy to lower frequencies as x/r increases. The effect of 
this increase will be discussed at the time the longitudinal 
scale of turbulence is considered. 


(a) Distance from jet. centerline y/r, 1.103. 

(b) Distance from jet centerline y/r, 1.00. 

Figure 24. — Comparison of longitudinal velocity correlations from 
one- and two-wire methods. 
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(a) to (f) Distance from jet centerline y/r, 1.00. 

(g) to (j) Distance from jet centerline y/r, 1.028. 

Figure 23. — Autocorrelations. Exit Mach number, 0.3; Reynolds number, 300,000. 


Spectral density, fi'i/), sec 


16 


REPORT 1292 — NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 



(a) Distance from jet centerline y/r, 1.206. 

(d) Distance from jet centerline y/r, 1.051. 

(g) Distance from jet centerline y/r, 0.846. 

Figure 25. — Spectral density curves. 


(b) Distance from jet centerline y/r, 1.154. 

(e) Distance from jet centerline ylr, 1.00. 

(h) Distance from jet centerline y/r, 0.794. 

Distance from jet nozzle x/r, 1.14; exit Mach 


(c) Distance from jet centerline y/r, 1.103 

(f) Distance from jet centerline y/r, 0.949. 

(i) Distance from jet centerline y/r, 0.743. 

number, 0.3; Reynolds number, 300,000. 


Spectral density, f(f) % sec 


INTENSITY, SCALF, AND SPECTRA OF TURBULENCE IN MIXING REGION OF FREE SUBSONIC JET 


17 



(a) Distance from jet centerline y/r, 1.00. 

(b) Distance from jet centerline y/r, 1.028. 

Figure 20. — Spectral density curves. Exit Mach number, 0.3; Reynolds number, 300,000. 
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(a) Distance from jet nozzle : c/r, 4.58. 

(b) Distance from jet nozzle x/r, 7.60. 

Figure 27. — Comparison of spectra obtained from a hot-wire anemom- 
eter and a microphone. Distance from jet centerline y/r, 1.00: exit 
Mach number, 0.3; Reynolds number, 300,000. 

There is a marked change in. the u spectra as the distance 
from the jet centerline y/r is changed. As the core of the 
jet is approached, a signal at approximately 8000 cps in- 
creases by more than an order of magnitude. This signal 
was believed to originate from a sound wave of that fre- 
quency. 

In figure 27 are compared microphone spectra and hot-wire 
spectra recorded for two different values of x/r with the 
microphone displaced in the ^-direction to remove it from 
the airstream. The presence of a strong signal at approxi- 
mately 8000 cps, which can be seen on both sets of records, 
shows that a sound wave is indeed the source of the extra- 
ne ous dist urba nc e . 

LONGITUDINAL SCALE OF TURBULENCE 
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(a) Various methods. 

(b) Various distances from jet centerline. 


Figure 28. — Variation of longitudinal scale with distance from jet 
nozzle. Distance from jet centerline y/r, 1.00: exit Mach number, 
0.3; Reynolds number, 300.000. 


curves. The results are shown in figure 28. A straight line 
has been fitted to the data by the method of least squares 
(fig. 28(a)). This straight line represents the variation of 
the longitudinal scale quite well for distances from the jet 
nozzle less than 16 radii. In figure 28(b) the results are 
plotted for downstream distances as great as 40 jet radii and 
show clearly the decrease in longitudinal scale as the distance 
from the jet nozzle is increased beyond 16 radii. 

The longitudinal scale of turbulence was evaluated by 
means of equation (13) and the spectral density curves of 
figures 25 and 26. It is obvious from equation (13) that 


( 21 ) 

Several schemes have been suggested to evaluate J( 0), the 
spectral density function for zero frequency. For example, 
von Karman (ref. 23) suggests that for small values of / an 
expression of the type 


The longitudinal scale of turbulence was determined 
by integrating the area under the longitudinal correlation 


7(f) 


h 

(l +»***/*)*'• 


( 22 ) 
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Figure 30. — Variation of longitudinal scale of turbulence across mixing 
zone. Exit Mach number, 0.3; Reynolds number, 300,000. 
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Distance from jet nozzle, x/ r 

Figure 31. — Variation of eddy size with distance from jet nozzle 

(eq. (23)). 


can be used to extrapolate to zero frequency. Corrsin and 
Uberoi (ref. 6) used this method to obtain the scale of tur- 
bulence. Figure 29 shows the result of fitting the data with 
a curve of the form of equation (22) for frequencies of 500 cps 
or less. It is seen that this method gives somewhat better 
results than equation (13). 

If this method of extrapolation is used, a value of the scale 
of turbulence can be obtained. It is emphasized, however, 
that the two values of the scale obtained (by fitting the corre- 
lation curve with the best exponential and by finding the best 
fit for the spectral density curve according to equation (22)) 
are not necessarily the same. This particular example shows 
the relative inadequacy of the method of equation (21) for 
obtaining the scale of turbulence. 

The method used in this report was simply to fair the best 
curve through the experimental data and estimate the maxi- 
mum value of J if). This value of J(j) * s called the J( 0) 
used in equation (21). This method gives values of £ x that 
agree well with the values obtained by integration of the 
eorrelograms, as is shown in figure 28(a), thus validating this 

method of selecting J( 0). The values of °Q X are plotted in 
figure 30. 

Examination of figures 17 and 28(a) shows that °Q X is 
somewhat more than twice as large as at x/r=S. In both 
figures, the values of and £ x for x/r= 1.14 fall below the 
least-squares lines. The scatter about that line is small ex- 
cept for this point. Inspection of figures 28(a) and (b) shows 
that it may be possible to fit the data with some combination 
of a parabola up to x/r= 16 and a straight line thereafter. 
This possibility suggests that the character of the turbulence 
changes at about that distance from the nozzle of the jet. 
In the mixing zone, the scale of turbulence is given by one 
relation, while the disappearance of the core results in an 
entirely different one. 

This fact is also shown by figure 30, which gives the relation 
of the longitudinal scale to the distance from the centerline 
of the jet. As long as there is a core, the relation is definitely 
one type, which changes markedly upon the disappearance 
of the core. This change, however marked it may be, does 
not influence the position of the largest eddies (£ Xtmax ) near 
the ylr value of 0.7 to 0.8 for values of x/r< 06. 

An analysis of the eorrelograms and spectrograms was 
made using the equation 

a _ _ x oO-~~P) (23) 

Az “ lnfl • ; 

The A scale of the turbulence was calculated by the methods 
outlined in appendix B. This A scale is different from that 
of the conventional scale £ x , but the over-all picture is much 
the same. Figure 31 shows the variation of this A scale 
with distance from the jet nozzle. Here the variation is 
linear with x/r , as was shown in figure 28(a). 


Spectral density, jr(/), sec 
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(a) Distance from jet centerline y/r, 0. 

(b) Distance from jet centerline y/r , 1.00. 

Figure 32. — -Spectral density curves for ^-component. Exit Mach number, 0.3; Reynolds number, 300,000. 
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(a) Distance from jet centerline y/r, 0. 

(b) Distance from jet centerline y/r, 1.00. 

Figure 33. — Velocity autocorrelations of u- and v-components. Exit Mach number, 0.3; Reynolds number, 300,000. 


CORRELATIONS AND SPECTRA OF u-COMPONENT OF 
TURBULENT VELOCITY 

Tiie correlations and spectra of the ^-components of tur- 
bulence were obtained in order to compare them with those 
of the ^-component. The spectra of the ^-component of 
turbulence at four points in the jet are shown in figure 32. 
The spectra of the ^-component of turbulence are different 
from the u - component in shape as well as in magnitude and 
frequency distribution. Figure 33 shows the autocorrela- 
tions of the v and the u velocity components at the same four 
points. These four curves indicate also the differences be- 
tween the u- and ^-components: 



r 


-► 



(b) 


r 


r 


If the values of „ and are related by the following 
equation (ref. 10) 


K 5 d7 ^ 

V.u-r 2 dx 


(24) 


then the turbulence is said to be isotropic. These u and v 
measurements, together with the v intensity measurements 
already discussed, show that the turbulence is not isotropic 
at distances from the jet less than 16 radii. 


(a) Band-pass width, 4200 t o 5800 (b) Band-pass width, 4800 to 5200 

cps; center f requeue} 7 , 5000 cps. cps; center frequency, 5000 cp.s. 
(c) Band-pass width, 900 to 1100 cps; center frequency, 100 cps. 

Figure 34. — Longitudinal velocity correlations. Distance from jet 
nozzle x/r, 1.14; distance from jet centerline y/r, 1.00; exit, Mach 
number, 0.3; Reynolds number, 300,000. 



Proportional to delay time, d/S , sec U Autocorrelation coefficient, 


22 


REPORT 1292 — NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 



( a ) No filter. (d) Band-pass width, 900 to 1100 cps. 

(b) Band-pass width, 4800 to 5200 cps. (e) Band-pass width, 480 to 520 cps. 

(c) Band-pass width, 2900 to 3100 cps. (f) Band-pass width, 190 to 210 cps. 

i ure 35. — Autocorrelations for specific band-pass frequencies. Distance from jet nozzle x/r, 4.58; distance from jet centerline y/r, 1.206; exit 

Mach number, 0.3; Reynolds number, 300,000. 



Figure 36. — Zero autocorrelation as function of frequency. Distance 
from jet nozzle x/r y 4.58; exit Mach number, 0.3; Reynolds number, 
300,000. 


CORRELATION WITH SPECIFIC BAND-PASS INPUTS 

As an additional test of the autocorrelation technique for 
obtaining longitudinal correlation coefficients, the hot-wire 
signals were passed through electronic filters set to pass only 
specific bands of frequencies before they were fed to the 
correlation computers. Figures 34(a) and (b) show the 
effect of varying the band-pass width. The correlations do 
not change materially as the width of the band pass is 
changed by a factor of 4, except that the narrow band pass 
produces a greater negative correlation (ref. 21). Figure 
34(c) shows the effect of changing the center frequency of 
the band pass. 

In figure 35 are shown the autocorrelations for the experi- 
ments where the hot-wire signals were passed through 
electronic filters before being fed to the correlation computer. 
This figure illustrates the effect of varying the center fre- 
quency of the band pass, as did figure 34(c), and also shows 
clearly how the presence of t lie filters changes the correla- 
tions as measured (e.g., cf. fig. 35(a), for which no filter 
was used, with the other parts of fig. 35, for which the filters 
were used). 
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(a) Distance from jet nozzle x/r, 1.14. (b) Distance from jet nozzle x/r, 2.2\). (c) Distance from jet nozzle x/r, 4.58. 

Figure 37. — Zero longitudinal correlation distance as function of frequency. Exit Mach number, 0.3; Reynolds number, 300,000. 


If the logarithm of the delay time for the first zero crossing 
of the filtered autocorrelation curve is plotted against the 
logarithm of the center frequency of the band pass (fig. 36), 
a straight line with a slope of — 1 results. If the logarithm 
of the wire separation for zero longitudinal correlation is 
plotted against the logarithm of the center frequency of the 
hand pass, a similar result should be obtained if the turbu- 
lent eddy has remained unchanged while passing the hot 
wires. That is, the eddy has been swept along unchanged 
with the velocity of the stream. Observation of figure 37 
shows that these curves are nearly straight lines with a 
slope of —1. This is additional evidence that the relation 
x=U i oc t holds very well for the turbulence patterns in the 
subsonic jet. Thus, the autocorrelation technique is a satis- 
factory one for studying the longitudinal correlations. 

VARIATION OF SCALE OF TURBULENCE WITH MACH AND/OR REYNOLDS 

NUMBER 

The effect of exit Mach and/or Reynolds number on the 
longitudinal scale was also investigated (fig. 38). As with 
the lateral scale (fig. 19), there was no variation of the 
longitudinal scale with Mach and/or Reynolds number. 



Reynolds number, Re 


Figure 38. —Effect of exit Mach and/or Reynolds number on longi- 
tudinal scale of turbulence. Distance from jet nozzle x/r, 4.58; 
distance from jet centerline y/r, 1.00. 

SUMMARY OF RESULTS 

Hot-wire-anemometer measurements of the turbulence 
parameters — intensity, scale, correlation, and spectra — in a 


subsonic jet have been reported for a range of exit Mach 
numbers from 0.2 to 0.7 and for Reynolds numbers from 
192,000 to 725,000, with the following results: 

1. Near the nozzle of the jet (distances less than 4 or 5 
jet diam downstream of the nozzle) the intensity of turbu- 
lence, expressed as a percent of core velocity, is a maximum 
at a distance of approximately 1 jet radius from the center- 
line and decreases slightly with increasing Mach and/or 
Reynolds number. At distances greater than 8 diameters, 
however, the maximum intensity moves out and decreases 
in magnitude until the turbulence-intensity profiles are quite 
fiat and approaching similarity. 

2. The lateral and longitudinal scales of turbulence are 
nearly independent of Mach and/or Reynolds number and, 
in the mixing zone near the jet, vary proportionally with 
distance from the jet nozzle. Farther downstream of the 
jet nozzle the longitudinal scale reaches a maximum and 
then decreases approximately linearly with distance. Near 
the nozzle the lateral scale is much smaller than the longi- 
tudinal and does not vary with distance from the centerline, 
while the longitudinal scale is a maximum at a distance from 
the centerline of about 0.7 to 0.8 of the jet radius. Farther 
downstream this maximum moves out from the centerline. 

3. An analysis of the correlograms and spectra, which 
differs from the ordinary statistical approach, yields a 
“scale” which, while different in magnitude from the con- 
ventional, varies similarly to the ordinary scale and is 
easier to evaluate. 

4. An an toco rrelator using a magnetic tape recorder and 
a special playback instrument was used to measure the 
velocity autocorrelations. These autocorrelations were con- 
verted to longitudinal correlations that agreed well with the 
directly measured longitudinal correlations. 

Lewis Flight Propulsion Laboratory 

National Advisory Committee for Aeronautics 
Cleveland, Ohio, April 24, 1956 
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APPENDIX A 


effective power band width 

displacement of heads of playback instrument 

spectral density function of u 2 

frequency 

calibration constant 

scale of turbulence 

exit Mach number 

number of db taken from strip chart 

correlation coefficient 

Reynolds number 

jet radius 

tape speed 

time 

delay time 

mean stream velocity 

fluctuating, components of velocity in x-, y-, and 
^-directions, respectively 

right-hand coordinate system with r-axis coincid- 
ing with jet centerline 

nondimensional parameter related to cut-off fre- 
quencies of spectrum by eq. (B5) 
fluctuating component of hot-wire voltage 
wire separation in x- and ^-directions, respectively 
characteristic eddy size 
angular frequency 


SYMBOLS 

Subscripts: 

b band of frequencies 

c core 

/ frequency 

1 lower cut-off 

loc local 

max maximum 

min minimum 

p peak 

ref reference 

t time 

tot total spectrum of all frequencies 

u upper cut-off 

w wire 

x longitudinal 

y lateral 

2 lateral 

1,2 points in stream 

Superscripts : 

average 

roo t -m ean -squa re 
* characteristic value 
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The analysis of the data from the hot-wire anemometers 
is most easily accomplished by means of the power spectrum 
of the turbulence and the autocorrelation technique, as 
explained in the EXPERIMENTAL PROCEDURE section. 
The details of this analysis are interesting also. The 
correlation fitting function 


e x/\x—p e 0*/A * 

^ 1 —p 

(14a) 

or 

1 ~P 

(14b) 

can be transformed to the spectral density function 
follows: 

Ju) as 

J(J)=tj f % cos j/ ch 

^ UlocJ 0 V loc 

(Bla) 

or 


*Rt cos 27r/£ d t 

(Bib) 

The resulting expression for J (j) is 


JU) 1+11 - (f) ■ 

4 ' ['X^Kf)] 

(17) 

This equation for y (J) has many interesting features. Con- 
sider, for instance, sketches (c) and (d). The typical cor- 
relogram shows a zero crossing (to) and a minimum value. 
These two quantities and the following equations for 

%,mm and o>* 


(B2) 

* 0 

(B3) 


give values for P and w* from the experimental correlogram. 
With these values the spectral density curve can be calcu- 
lated from equation (17). 

An alternate procedure is to work with the experimental 
spectral density curve. If the low-frequency part of the 
curve has been obtained with sufficient accuracy, the values 
of the low- and high-frequency cut-offs (f t and /J can be 
obtained from the spectrum. These two values with the 
equations 


fi_ P 

ju ( 1 + 0) 2 


(B4) 


and 


_(l+0)«* 

Ju 2tt 


(BS) 


give values for p and co* from the experimental spectrum 
With these values the correlogram can be calculated from 
equation (14a) or (14b). 

In actual practice the low-frequency part of the spectrum 
is seldom defined with sufficient accuracy to use the alternate 
method proposed here. But the eorrelograms are generally 
well defined and are easily used as outlined. 

If a family of curves of f (j) for a series of values of p is 
plotted as shown in figure 39, spectrum fitting curves result 
that differ from experimental curves only by multiplicative 
constants. This family of curves can be used as a grid and 
placed over the experimental data and a best fit obtained. 
The value of 0, therefore, is easily determined. With this 
value of P and the zero crossing of the correlogram t 0} the 
value of co* can be obtained from equation (B3). 

With the values of p and co* the spectrum can be calcu- 
lated by means of equation (17), or A x can be evaluated by 
means of equation (16). 



(c) 



(d) 
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Figure 39. — Spectrum fitting curves (eq. (B2)). 
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